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Abstract 

Mortar methods have recently been shown to be well suited for isoge¬ 
ometric analysis. We review the recent mathematical analysis and then 
investigate the variational crime introduced by quadrature formulas for 
the coupling integrals. Motivated by finite element observations, we con¬ 
sider a quadrature rule purely based on the slave mesh as well as a method 
using quadrature rules based on the slave mesh and on the master mesh, 
resulting in a non-symmetric saddle point problem. While in the first 
case reduced convergence rates can be observed, in the second case the 
influence of the variational crime is less significant. 


1 Introduction 

Isogeometric analysis, introduced in 2005 by Hughes et al. in [IHCB051 . is a fam¬ 
ily of methods that use B-splines and non-uniform rational B-splines (NURBS) 
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as basis functions to construct numerical approximations of partial differential 
equations (PDEs). With isogeometric methods, the computational domain is 
generally split into patches. Within this framework, techniques to couple the 
numerical solution on different patches are required. To retain the flexibility of 
the meshes at the interfaces, mortar methods are very attractive. 

Mortar methods are a popular tool for the coupling of non-matching meshes, 
originally introduced for spectral and finite element methods [BMP941IBM971 
IBen99j . They were successfully applied in the context of isogeometric analy¬ 
sis |HB12I IASWB141IDVK14] . A mathematical analysis enlightening the use 
of different dual spaces was recently presented in [BBWWT5] , In this paper, 
starting from these latter results, we focus on one particular challenge in the 
realization of a mortar method, namely, the evaluation of the interface integrals 
which contain a product of functions defined on non-matching meshes. 

This article is structured as follows. In Sec.[2j we recall the basics of isogeo¬ 
metric mortar methods. In Sec. [3] we consider a review of numerical quadrature 
for mortar integrals such as additional aspects in the case of isogeometric anal¬ 
ysis, illustrated by numerical results in Sec.[4j 

2 Isogeometric Mortar Methods 

In the following, we briefly present isogeometric mortar methods, for more de¬ 
tails we refer to |BBWW15] . After stating the problem setting, we review 
isogeometric parametrizations, describe the domain decomposition into several 
NURBS patches and finally discuss suitable coupling spaces. 

Let fl C M d , d the dimension being 2 or 3, be a bounded domain, a, /3 £ 
a > «o > 0 and /3 > 0. We consider the following second order elliptic 
boundary value problem with homogeneous Dirichlet conditions 


div(aVw) + j3u = f in 12, 

u = 0 on dflo = 912. 


(la) 

(lb) 


We assume a and /? to be piecewise sufficiently smooth. 

2.1 Isogeometric Parametrization 

Here, we present isogeometric concepts and no tations used throughout the pa¬ 
per, and refer to the classical literature [PT97I lBBdVC + 06l ISch07l ICHB09] for 
more details. 

Let us denote by p the degree of the univariate B-splines and by S an open 
univariate knot vector, where the first and last entries are repeated (p+1)-times, 
i.e., 

^ = {0 = £l = • • ■ = £p+l < £p+2 fs • • • fs £ri < £n+l = ■ • • = £n+p+l = 1}- 

Let us define Z = {£ 1 , £ 2 , • • •, Ce} as the knot vector without any repetition, 
also called breakpoint vector. For each breakpoint Q of Z, we define its mul¬ 
tiplicity rrij as its number of repetitions in S. The Cox-de Boor algorithm, 
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see [CHB09) , defines n univariate B-splines Bf((),i = 1 ,n, based on the uni¬ 
variate knot vector S and the degree p. We denote by S p ( 3) = span{13f(£), i = 
1, ..., n} the corresponding spline space. The smoothness of B-splines is defined 
by the breakpoint multiplicity. More precisely, each basis function is C p ~ m:i at 
each (j e Z. 

To define multivariate spline spaces, we introduce the multivariate knot vec¬ 
tor S = (Si x S 2 x ... x 3,j), and for simplicity of notations assume in the 
following that the degree is the same in all parametric directions and denote it 
by p. Multivariate B-splines B p (£,) are defined by tensor product of univariate 
B-splines for each multi-index i 6 I = {(*i, ..., id) : 1 < is < ng}. We denote 
by S p ( S) the corresponding spline space in the parametric domain. 

Given a set of positive weights we define NURBS functions N ?(£) as 
rational functions of B-splines and the weight function W = Xaei w i (C), and 
set N p (a) as the multivariate NURBS space in the parametric domain. 

For a set of control points C, £ M. d , i £ I, we define a parametrization of a 
NURBS surface (d = 2) or solid (d = 3) as a linear combination of NURBS and 
control points 

F(C) = ^QfVf(C), 

iei 

and assume the regularity stated in |BDVBSV141 , Assumption 3.1]. 

The knot vector S forms a mesh in the parametric space fl. We define the 
physical mesh A4 as the image of this parametric mesh through F, and denote by 
O its elements. The /i-refinement procedure, see [BDVBSV141 Section 2.1.3], 
yields a family of meshes denoted Mh, each mesh being a refinement of the 
initial one. We assume quasi-uniformity for each mesh. 


2.2 Description of the Computational Domain 


Let the domain fi be decomposed into K non-overlapping domains fh, i.e., 


K 

fl = Ofc, and U,; D ft,- = 
k= 1 


U l 3 - 


Each subdomain is a NURBS geometry, i.e., there exists a NURBS parametriza¬ 
tion Ffc based on a knot vector S* and a degree pk, see Sec. |2.1| such that fh 
is the image of the parametric space U = (0, \) d by F k . 

For 1 < fci, < K , fci fc 2 , we define the interface as the interior of the 
intersection of the boundaries, i.e., 7 klko = dilk 1 D d£lk 2 , where 7 k x k 2 is open. 
Let the non-empty interfaces be enumerated by 7 ;, l = 1, ..., L, and let us 
define the skeleton T = U^Li li as the union of all interfaces. For each interface, 
one of the adjacent subdomains is chosen as the master side and one as the slave 
side. This choice is arbitrary but fixed. We denote the index of the former by 
m(Z), the index of the latter one by s(l), and thus 7 y l = 2 m p) H dCl s (iy On 
the interface 7 ;, we define the outward normal n/ of the master side and 

du 

denote by —— the normal derivative on 7 ; from the master side. 
oni 
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Figure 1: Geometrical conforming case (left) and slave conforming case (right). 


We assume that for each interface the pull-back with respect to the slave 
domain is a whole face of the unit d-cube in the parametric space, which we call 
a slave conforming situation, see the right setting in Fig. |T] If we also assume 
that the pull-back with respect to the master domain is a whole face of the unit 
d-cube, we are in a fully geometrically conforming situation, see the left picture 
of Fig. [T] 

For each Q k , we introduce the space Hl(fl k ) = {v k £ H 1 (n k ),v k \ gflngn = 
0}, where we use standard Sobolev spaces, as defined in [Grill) , endowed with 
their usual norms. For any interface 7 ; C <9f2 s p), we define by H 0 ' Q ( 7 ;) C 
id 1 / 2 (9n s (/)) the space of all functions that can be trivially extended on <9fl s p) \ 
7 ; by zero to an element of id 1 / 2 (9fl s (p). Note that id 1 / 2 (9n s (p) is the trace 

space of Jd 1 (f2 s (/)). The dual space of H^{^i) is denoted it _ 1 / 2 ( 7 ^). In order to 
set a global functional framework on Q, we consider the broken Sobolev spaces 
V = Ylk=i H *(^k), endowed with the broken norm ||u||^ = J2k=i IMIaRa*)’ 
and M = J]f = i # _1/2 ( 7 *)- 

The mortar method is based on a weak coupling between different subdo¬ 
mains. Each subdomain is discretized independently and a weak coupling is 
performed on each interface. From now on, we assume that jumps of a and /? 
are solely located at the skeleton, and we define the linear and bilinear forms 
a: V x V —> I and /: V —> R, such that 


K f 

a(u,v) = E / aS7u ■'Vv + f3u v dx, 

fc=i 



fv dx. 


We remark that the standard weak formulation of 0 . where no weak coupling 
is necessary and which is uniquely solvable, reads as follows: Find u £ i?o(fl), 
such that 


a(u, v ) = f(v), v £ Hq(Q). ( 2 ) 

2.3 Isogeometric Mortar Discretization 

In the following, we set our non-conforming approximation framework. On 
each subdomain £l k , based on the NURBS parametrization, we introduce the 
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approximation space 14 ,h = {i’k = Vk £ N Pk ( H^)}. We recall that 

under the assumptions on the mesh A ik,h and on the parametrization F^, this 
NURBS space has optimal approximation properties, see, e.g 


BBdVC+06 


On LI, we define the discrete product space 14 = JlfcLi 14 ,h C V, which forms 
a H 1 (Q) non-conforming space discontinuous over the interfaces. We denote in 
the following the maximal mesh size h = max^ hk as the mesh parameter. 

On the skeleton T, we define the discrete product Lagrange multiplier space 
M h as M h = nf = .j M/j,, where denotes one of the two following choices 
which were shown to be well suited in |BBWWl5| . The first choice M® h is the 
spline space of degree p s (i), defined on the interface 7 1 based on the interface 
knot vector of the slave body Sl a a). Note that in the presence of any cross point, 
a suitable modification, e.g., a local degree reduction as presented in [BBWW151 , 
Section 4.3], has to be applied. The alternative choice Mf h is an order (p s (p — 2) 
spline space defined on the interface 7 1 based on the interface knot vector of the 
slave body U s p) for which the definition requires the trace space of V s ^ t h to be 
a subset of C 1 ^). Moreover, a third dual space M / , which is a spline space of 
degree (p s (p — 1), could be considered. This choice is not inf-sup stable, so we 
do not consider it any further. 

The saddle point formulation of the isogeometric mortar method reads as 
follows: Find (uh,Xh) £ 14 x Mh, such that 


a(u h ,v h )+ b(v h ,X h ) = f(v h ), v h &V h , 
b(u h ,p, h ) = 0, Hh&M h , 


(3a) 

(3b) 


where b(v,i Li) = / 7i mMz d * 7 an d [•]/ denotes the jump from the master 

to the slave side over 7 /. We note that the Lagrange multiplier A h gives an 
approximation of the normal flux across the skeleton. 

It is well known from the theory of mixed and mortar methods that two 
abstract requirements on each interface guarantee the method to be well-posed 
and of optimal order, see |Ben99] . Namely, an appropriate approximation order 
of the dual space and a uniform inf-sup stability between the primal space and 
the dual space. Note that for simplicity of notations, we assume the same type 
of dual space to be used for all interfaces. The following theorem is shown 
in |BBWW15] and guarantees a-priori bounds. 

Theorem 1 Let 9 = 0 if t and 8 = 1/2 if Mf, i = ; . For 

u £ H a+l {FL), 1/2 < tr < min k,i{Pk ~ 0), the solution of ([2]) and {uh,\h) the 
non-conforming approximation, see it holds 


K 


j^\\ u ~ M ^lll 2 ( 0 ) + ll u ~ u h\\v ^ 4 1 ^2 hl a \\ u \\'k^+ 1 (n k ) 


k -1 


for the primal solution and 


L n k 

J2w " 2 


1=1 


a dn, 

1 k=1 


r M 2 H « +H n k) 
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for the dual solution. With 0 < C < oo a generic constant that is independent 
of the mesh sizes but possibly depends on pk- 

We highlight that while for both pairings the inf-sup stability is satisfied, the 
approximation order of the lower order dual space Mf h is only close to optimal. 
Indeed, in this case an order of \J~h is lost in the convergence order. 

3 Mortar Integrals 

To evaluate the bilinear form 6 (u,/i), we need to evaluate for each interface 7 ; 
the mortar integrals p, v + da and f ^tu~der, where v + denotes the trace of 
v from the master domain fl m (;) and v~ the trace of v from the slave domain 
^ s(l )• To simplify the notation, let us restrict ourselves to the case of one single 
interface and drop the index l in the following. 

One particular challenge in the realization of a mortar method is the evalua¬ 
tion of the first interface integral due to the product p, v + of functions which are 
defined on non-matching meshes, see |BF04] for a method to bypass it in a finite 
element/wavelet context. Any quadrature rule based on the slave mesh does not 
respect the mesh lines of the master mesh and vice versa for a quadrature based 
on the master mesh. 

It is obvious that the use of a suitable quadrature rule based on a merged 
mesh, i.e., a mesh which respects the reduced smoothness of the master and 
slave functions at their respective mesh lines, leads to an exact evaluation of the 
integral. However, the construction of this auxiliary mesh commonly named 
segmentation process is challenging, especially in the three dimensional case 
since the shape of the elements is not unique and difficult to determine, see, 
e.g., |ML00l IPL041 lPus041 IHB121 IDFTH14| . Note that in an isogeometric con¬ 
text the merged mesh needs to be constructed in the physical space and then 
pulled back to the parametric space for each subdomain. The complexity of 
constructing such a mesh becomes even more severe in the case of non-linear 
and time-dependent problems, where the relative position of the meshes changes 
in every time or load step which implies to recompute the merged mesh at every 
step. 

Due to this computational complexity, it has been seen very appealing to 
use a higher order quadrature rule either based on the slave mesh or on the 
master mesh, see [;FW051 ITFW09IIDLTWZTT] for some applications in finite 
element and isogeometric analysis contexts. However in the finite element case, 
early results in [CLM971 IMRW02] showed that this strategy does not yield 
optimal methods. More precisely, in the case the master mesh is chosen, the 
best approximation error is affected, while in contrast in the case the slave 
mesh is chosen it is the consistency error. Numerical results confirmed the 
lack of optimality with the master integration approach, while with the slave 
integration approach reasonable results were obtained although not optimal in 
terms of the Lagrange multiplier norm. 

Due to the global smoothness of splines, one could expect the sensitivity with 
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respect to the quadrature rules for isogeometric methods to be less than for finite 
element methods. In the mortar context, according to the finite element results, 
it seems interesting to consider a slave integration rule. And, in case of maximal 
regularity, i.e., 14 ,h C C p ' c_1 (f4) one also might expect the quadrature error 
on a non-matching mesh to be significantly smaller than in the finite element 
case. 

Let us denote the quadrature rule based on the boundary mesh of the slave 
domain as ]P_, be., f 7 l^v + da « '}2-l jlv+ - We precise that in the examples a 
Gaussian quadrature rule with a various number of points is used. The mortar 
method with pure slave integration is obtained by evaluating all interface inte¬ 
grals in a using this quadrature rule, i.e., the discrete system reads as follows: 
Find (uh, A/i) £ 14 x Mh, such that 

a{u h ,v h ) + ~ v h)^h = f{vh), v h £ 14 , 

y: (ui -u^)nh=o, ^ € m h . 

The notation ~ is used to stress the difference to the discrete solution with exact 
integration. 

In the next section, we present numerical examples which show severe distur¬ 
bances even in the isogeometric case. Hence, even though the global smoothness 
of the integrated function is increased compared to the finite element case, a 
non-matching integration approach reduces the convergence order drastically. 

Moreover, we consider an alternative approach which was proposed in [CLM971 
IMRW02~| using both integration rules. Additionally denoting a quadrature 
rule based on the boundary mesh of the master domain f l m , this approach 
which results in a non-symmetric saddle point problem, reads as follows: Find 
(u h , X h ) £ 14 x M h , such that 

a(u h ,v h ) + y; Vh^h -^2_ v h^h = f(vh), Vh £ 14, 

- u^)fi h = 0, Uh £ M h . 

The non-symmetric saddle point problem, which corresponds to a Petrov-Galerkin 
approach in the primal formulation, was motivated by different requirements 
for the integration of the primal and dual test functions. Numerical examples 
showed error values very close to the case of exact integration, but we note that 
from the theoretical side even the well-posedness of the non-symmetric saddle 
point problem remains unclear. In the next section, we present numerical ex¬ 
amples which show that also in an isogeometric context the results are generally 
close to those from the exact integration case. 

4 Numerical Results 

In this section, we consider two-dimensional and three-dimensional settings in 
order to observe the effects of inexact quadrature rules, as presented in Sec. [3] 
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Figure 2: Left: Primal solution on Q. Right: Lagrange multiplier along the 
interface. 


on the optimality of the mortar method. We first set the problem settings, and 
then give the results of several studies. 

4.1 Two-dimensional Example 

As a first example, let us consider the Poisson problem —A u = f solved on 
the domain = ( 0 , 1 ) x (— 1 , 1 ) which is decomposed into two patches by the 
interface 7 = {(a:, y) £ O, y = 0}. The upper domain is set as the slave domain. 
The internal load and the boundary conditions are manufactured to have the 
analytical solution 



The normal derivative on the interface is given by du/d n(*) = 27rcos ( 7 ra;), see 
Fig.|] Neumann conditions are applied on the left and right boundary parts, 
such that no cross point modification is necessary. 

Regarding the meshes, we consider three different cases, presented in Fig. [3] 
In the first two cases, the initial master mesh is a refinement of the initial slave 
mesh. The initial slave mesh consists of just one element. In the case Ml, one 
uniform refinement step is applied to build the master mesh, in the case M2 two 
uniform refinement steps. Case M 3 was chosen such that at no refinement level 
parts of the slave and master boundary meshes do coincide. The initial interior 
knots of the slave domain were chosen as { 7 r/ 10 ,1 — 7 t/T} in both parametric 
directions, yielding 9 elements. The initial master mesh consists of four uniform 
elements. 

In the following, we provide different numerical error studies. We note that 
the inter-element smoothness of the dual functions can influence the accuracy 
of the quadrature based on the master mesh, but not the one based on the slave 
mesh. Therefore for the slave integration approach, the equal order pairing with 
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Figure 3: Different meshes at mesh refinement level 1. From the left to the 
right: Ml to M3. 


maximal smoothness is considered, i.e., Mh = M° C C p ~ 1 ( 7 ), while for the non- 
symmetric approach we vary the dual degree. In all cases, the primal L 2 (fl) and 
the dual L 2 ( 7 ) errors are computed by a comparison with the analytical solution 
stated above. 

4.1.1 Slave Integration Approach 

Firstly, we consider the case M3, see Fig. [3] to measure the impact of the 
integration error in a general situation. A numerical error study is provided 
in Fig. [4] for a different number of Gauss points and different spline degrees. 
It can clearly be seen that the primal and dual solutions are both affected 
by the inexact quadrature, leading to non-optimal methods. In all cases, the 
same characteristic behavior can be seen. Up to a certain refinement level, the 
results with inexact quadrature rules coincide with the ones with no quadrature 
error. Then, at a certain refinement level, the convergence order is reduced and 
the error is significantly larger than the exact integration one. The starting 
disturbance threshold is different for the primal and dual solutions as well as it 
differs for a different amount of quadrature points. Moreover, in this situation 
the higher order splines are more disturbed by the numerical quadrature than 
the lower order splines. 

In almost all cases of Fig. [4] we observe poor approximation results and a 
reduced convergence order which is numerically independent of the spline degree. 
Especially, the rate of the L 2 ( 7 ) dual error is very low, so let us consider the final 
numerical convergence rate in more details. In Table [T] estimated convergence 
orders for degree 5 and cases Ml and M2 are given. We notice that the dual 
L 2 ( 7 ) rate breaks down to an order of 1/2, while the L 2 (Q) primal rate lies 
about 3/2. 

Secondly, we consider an even more simple situation to show that even then 
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Figure 4: 2D results - L 2 primal (left) and dual (right) error curves for the case 
M3: equal order pairings with p = 1,3,5 (from top to bottom) for the slave 
integration approach and different quadrature rule orders. 
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primal error 

dual 

error 

case Ml 

case M2 

case Ml 

case M2 

quad, rule order 0 

1.63 

1.74 

0.50 

0.50 

quad, rule order 1 

1.63 

1.54 

0.50 

0.50 

quad, rule order 2 

1.63 

1.55 

0.50 

0.50 

quad, rule order 3 

1.63 

1.58 

0.50 

0.50 

quad, rule order 4 

1.63 

1.56 

0.50 

0.50 

quad, rule order 5 

1.63 

1.50 

0.50 

0.50 


Table 1: 2D results - Last estimated order of convergence of the primal and dual 
L 2 errors for the cases Ml and M2: pairing P 5 — P 5 for the slave integration 
approach and different quadrature rule orders. 


the impact of the slave integration is noticeable. Let us focus on the cases Ml 
and M2, see Fig. [3] for which the master mesh is a refinement of the slave 
mesh. See Fig. [5] for a comparison of results between the cases Ml and M2 for 
a spline degree p = 3. We note that the low convergence order of the primal 
and dual solutions, as remarked above, already appear in this simple context. 
Moreover, for a fixed number of slave elements, the error is increasing with the 
number of master elements. This is expected as there are more points of reduced 
smoothness which are not taken account by the quadrature rule. 

Thirdly, we have additionally compared the case M 3 with a similar situation 
in which the master and slave roles are inverted. The results also show that the 
integration error is increasing with the increase of the master element number. 
Thus, in accordance to the practical applications, in a slave integration context 
it seems worthwhile to choose the slave domain as the finest one. 

Moreover, it can be observed that on coarse meshes using the slave integra¬ 
tion method it is possible to recover the accuracy of the optimal mortar method 
simply by increasing the number of quadrature points, see Fig. [g] However, it 
has also been shown that the number of necessary quadrature points is drasti¬ 
cally increasing with the refinement level. It can easily be seen that the number 
of Gauss points gets soon impracticably large, see the right picture of Fig. [6] 
Furthermore, in several cases, the disturbance to the mortar method has been 
observed to be more severe for higher order functions. 

4.1.2 Non-symmetric Approach 

The non-symmetric saddle point problem based on the two different quadrature 
rules, see Sec. [3] was introduced to overcome the non-optimality of the pure 
slave integration approach in a finite element context. Due to the suboptimal 
results seen in the previous section, it is also interesting to consider it in an 
isogeometric context. 

Firstly, we consider same degree pairings. In almost all tested cases the 
results of the non-symmetric approach are comparable to the results of the 
exact integration case. However, we note that differences could still be seen 
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Figure 5: 2D results - L 2 primal (top) and dual (bottom) error curves for the 
cases Ml (left) and M2 (right): pairing P 3 — P 3 for the slave integration 
approach and different quadrature rule orders. 




Figure 6: 2D results - L 2 primal (left) and dual (right) errors at refinement level 
number 6 as a function of the quadrature rule order for the case M3: pairing 
P 3 — P 3 for the slave integration approach. 
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Figure 7: 2D results - L 2 primal (left) and dual (right) error curves for the case 
M3: equal order pairing p = 1 for the non-symmetric approach and different 
quadrature rule orders. 


in some cases. For example, for a degree p = 1 in the case M3, we obtained 
a non-optimal method, see in Fig. [7] the corresponding primal and dual error 
curves. Note that we do not show any curves in the cases where no disturbance 
is observed. For example for degree p = 5 we observed convergence almost up 
to machine precision without any remarkable difference compared to the exact 
integration case. 

Secondly, we consider dual spaces with lower degrees than the primal ones. 
Note that in (BBWWI5] stability for these pairings was only observed if the 
primal and the dual degrees have the same parity. Similar to the equal order 
case, the dual error is not affected. In Fig. [ 8 ] primal error curves are shown 
for all stable different degree pairings up to a primal degree p = 4. We note 
that theoretically, we expect sub-optimal primal error rates even in the exact 
integration case, although often improved convergence rates were observed. For 
a dual degree p—2k, k £ N, we can expect a convergence of order Oih p ~ 2k+b / 2 ) 
in the L 2 (fl) norm, see the dashed lines in Fig. [ 8 ] For the P4 — P2 and P3 — PI 
pairings, we observe slight differences compared to the exact integration results, 
but note that the convergence rate is not far from the theoretical expectation. 
The situation is different for the P4 — P0 and P2 — P0 pairings, for which the 
rate is more disturbed and even below the theoretical expectation. This could 
be explained by the discontinuity of the dual basis functions, which introduces 
large errors in the integration approximation done with a rule based on the 
master mesh, which does not respect these discontinuities. 

4.2 Three-dimensional Example 

As a second example, we consider a three-dimensional problem with a curved 
interface. Precisely, we consider the Poisson problem —A u = f on the do¬ 
main = ( 0 , l) 3 , which is divided into two patches by the interface 7 = 
{(x,y,p(x)),(x,y) G (0,1) 2 }, with p{x,y) = 1/8 (1 + x)(l + y 2 ) + 1/5, see 
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Refinement level number Refinement level number 




Refinement level number Refinement level number 


Figure 8: 2D results - L 2 primal error curves for the case M3: different order 
pairings for the non-symmetric approach and different quadrature rule orders. 
Top left: P2 — PO. Top right: P3 — PI. Bottom left: PA — P2. Bottom right: 
PA - PO. 
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Figure 9: Meshes at refinement level 1 (left) and the slave domain (right) illus¬ 
trating the curved interface. 

Fig.m The bottom domain is set as the slave domain. The internal load and 
the boundary conditions are manufactured to have for analytical solution 

u(x,y,z ) = cos(2TTx)cos(2ny)sin(2TTz). 

Note that due to the curved interface, the normal derivative has a complicated 
form, but is still explicitly computable. Neumann conditions are applied such 
that no cross point modification is necessary. The initial master mesh has 
8 uniform elements, while the initial slave mesh has 8 elements given by the 
breakpoint vector {0,7r/5,1} in each direction. In the following, we provide 
some numerical error studies, considering the slave integration approach as well 
as the non-symmetric approach. 

The obtained results are in accordance with the two-dimensional results for 
both approaches. In Fig. [lOj the disturbance for the slave integration approach 
is shown for the P4 — PA pairing. Although not shown here, we note that the 
results for the P2 — P2 and P3 — P 3 pairing have a similar behavior. The 
non-symmetric approach does not lead to reduced rates considering equal order 
pairings, i.e., M^ = M®, on the refinement levels we considered. As previously, 
using a lower order dual space, a difference to the exact integration case can be 
seen. See Fig. [Tljfor the disturbance in the primal variable of the P 3 — PI and 
P4 — P2 pairings. 


5 Conclusion 

In this article, after reviewing optimal isogeometric mortar methods, a study 
on the possibility to approximate the mortar integrals by efficient numerical 
quadrature rules was performed. 
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Figure 10: 3D results - L 2 primal (left) and dual (right) error curves for the 
pairing P4 — P4, for the slave integration approach and different quadrature 
rule orders. 




Figure 11: 3D results - L 2 primal error curves for the pairings P3 — PI (left) 
and P4 — P2 (right), for the non-symmetric approach and different quadrature 
rule orders. 
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To exactly integrate the product of functions defined on non-matching meshes, 
as in the mortar integrals, it is necessary to construct a merged mesh. Since this 
construction is of a high complexity, it would be desirable to use a quadrature 
rule based on the slave mesh only. However, numerical examples show a sig¬ 
nificant disturbance to the mortar method, especially for higher order splines. 
Especially the convergence rate of the Lagrange multiplier is reduced to 1/2 
and less. While the method improves by increasing the number of quadrature 
points, the amount of points necessary to obtain nearly optimal results is not 
predictable. 

To overcome these difficulties, we have considered a non-symmetric saddle 
point problem based on both master and slave integration rules, which was previ¬ 
ously introduced in the finite element context. Numerical examples demonstrate 
the possibility to reach the accuracy given by an exact integration strategy, al¬ 
though this it is not ensured for all cases. 
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